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We revisit the phenomenon of spinodals in the presence of quenched disorder and develop a 
complete theory for it. We focus on the spinodal of an Ising model in a quenched random field 
(RFIM), which has applications in many areas from materials to social science. By working at zero 
temperature in the quasi-statically driven RFIM, thermal fluctuations are eliminated and one can 
give a rigorous content to the notion of spinodal. We show that the latter is due to the depinning 
and the subsequent expansion of rare droplets. We work out the associated critical behavior, which, 
in any finite dimension, is very different from the mean-field one: the characteristic length diverges 
exponentially and the thermodynamic quantities display very mild nonanalyticities much like in a 
Griffith phenomenon. From the recently established connection between the spinodal of the RFIM 
and glassy dynamics, our results also allow us to conclusively assess the physical content and the 
status of the dynamical transition predicted by the mean-field theory of glass-forming liquids. 

PACS numbers: 


The instability or failure of a metastable state in the 
presence of quenched disorder plays a key role in a 
wide spectrum of out-of-equilibrium situations, where it 
leads to extreme events such as macroscopic avalanches, 
shocks, ruptures or crises. Such “spinodal transitions” 
can be found for instance in the hysteretic behavior of 
disordered magnets [l|, 0 , capilla^ condensation in dis¬ 
ordered mesoporous materials |3|,1^, failure in amorphous 
materials [^0, and in a variety of social and economic 
phenomena [7[. Whereas spinodal transitions are well 
understood in “clean” systems, where the limit of stabil¬ 
ity of a thermodynamically metastable state is associated 
with the emergence of “soft modes” , quenched dis¬ 
order introduces new phenomena and drastically changes 
the physical behaviour: it pins the soft modes so that the 
physics is dominated by avalanches and bursts in place 
of long wave-length fluctuations. 

In this work we present a complete theory of spinodal 
transitions in the presence of quenched disorder. In ad¬ 
dition to its wide relevance from materials to social sci¬ 
ence, one main motivation to focus on this problem is 
its connection with glassy dynamics 11|, as unveiled by 
recent advances in the theory of glass formation. Indeed, 
in what now appears to be the proper mean-field (MF) 
theory of the glass transition [1^ [l3| , the liquid as it is 
cooled first undergoes a dynamical transition, which is 
essentially that described by the mode-coupling theory 
(MGT) This transition is found to be in the same 

universality class as the spinodal of the random-field Ising 
model (RFIM) [HI- The fact that the dynamical tran¬ 
sition is the point where metastable states emerge (or 
loose their stability if one comes from low temperature) 
and that these states correspond to disordered particle 
arrangements provides the physical basis for this result. 


From field-theoretical arguments, the upper critical di¬ 
mension above which the critical exponents take their 
MF value, i.e., here, the MGT predictions, was argued 
to be due = 8 iH 3, as indeed also proposed for dis¬ 
ordered spinodals [l^ ■ The theoretical understanding of 
the spinodal of the RFIM in finite dimensions is therefore 
the crucial missing step to conclusively assess the physi¬ 
cal content of the MGT and the status of the dynamical 
transition predicted by the MF theory of glass-forming 
liquids. 

It is well known however that rigorously speaking a spin¬ 
odal cannot in general exist in finite dimensions, even 
above the upper critical dimension. This is due to the 
presence of thermal fluctuations that destroy metasta¬ 
bility through nucleation phenomena. What could at 
best remain are therefore vestiges of the singularity in the 
form of a crossover behavior (as, e.g., seen between spin¬ 
odal decomposition and nucleation-and-growth in demix¬ 
ing transitions). However, in several cases, and for the 
RFIM, the fluctuations induced by the quenched disor¬ 
der play a central role. The long-distance physics can 
then be directly studied at zero temperature [17|. This 
in particular implies that a spinodal can now be rigor¬ 
ously defined beyond MF. 

In the following, we focus on the RFIM and present a 
complete theory of its zero-temperature spinodal in all 
dimensions. The Hamiltonian of the model reads 


n = -jY^ 


( 1 ) 


<v> 


where Ising spins S'i = ±1 are placed on the vertices of 
a d-dimensional hypercubic lattice, < ij > denotes dis¬ 
tinct nearest-neighbor pairs, the hi’s are random helds 
independently taken from a Gaussian distribution of zero 
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mean and variance = R^, and H is the external mag¬ 
netic field. The coupling J is taken as the unit of energy 
and temperature (J, fcs = 1). 

Studying the zero-temperature spinodal of the RFIM 
means following its hysteretic, metastable behavior: one 
changes H quasi-statically and let the system evolve un¬ 
der a r = 0 Metropolis dynamics: a spin is flipped as 
soon as the local field at its site, J2j\i + hi + H, 
changes sign, whatever the location of this spin [2, [l^ 
[see also the supplemental material (SM) [39|]. There 
is a critical value of the disorder R = Rc that sepa¬ 
rates a high-disorder regime where the hysteresis loop 
in the magnetization-versus-field plane is continuous and 
a low-disorder regime where a macroscopic magnetiza¬ 
tion jump is observed 0. In the former, the avalanches, 
which are discontinuous changes of the system’s config¬ 
urations and are typical of the evolution of disordered 
systems at T = 0, are of limited size whereas in the lat¬ 
ter an infinite, system-spanning and compact, avalanche 
is found. For R < Rc, the value of the applied field (the 
“coercive field”) at which the system spanning avalanche 
takes place can then be considered as the spinodal tran¬ 
sition of the model: It is the limit of stability of the pos¬ 
itively magnetized phase when the field is decreased or 
of the negatively magnetized one if the field is increased. 
(For R > Rc, there is no transition at all.) 

The key questions we address are: Under what condi¬ 
tions, if any, does this well-defined RFIM spinodal, which 
is in the same universality class as the MCT dynamical 
transition of glass-forming liquids, display critical behav¬ 
ior? What is the physical mechanism governing the spin¬ 
odal transition and its long-distance universal behavior? 
A naive response, supported by some numerical work 
0 , is that the spinodal should be MF-like and critical, 
i.e., with diverging correlation length and susceptibility, 
above an upper critical dimension due = 8 and noncriti- 
cal, i.e., with finite correlation length and susceptibility, 
below. By combining analytical and numerical investiga¬ 
tions, we show that this is not what happens. The critical 
behavior predicted by MF theory is absent in any finite 
dimension, even above due = 8. It is destroyed not only 
by thermal fluctuations, which was already well known, 
but also by rare athermal nonperturbative fluctuations. 
There is still a transition at T = 0, but it greatly differs 
from its MF counterpart. 

We start with a surprising outcome of computer sim¬ 
ulations for 3-dimensional cubic lattices and periodic 
boundary conditions (PBC). For R < Rc ~ 2.2 a spin¬ 
odal takes place, but the value of the coercive field, 
Hc{R), increases as one decreases the disorder below the 
critical value Rc'. see for instance Fig. 1 of Ref. [l^ . 
Figs. 8 and 16 of Ref. [l^, Fig. 1 of Ref. 211, and 
the triangles in Fig. [TJi (obtained for linear size L = 90 
and averaging over 2400 samples). This is counterin¬ 
tuitive as one expects the opposite trend, with weaker 
disorder in general reducing the width of the hystere¬ 


sis loop 22|. This is even more puzzling if one studies 
the R ^ 0 limit. In this case one can show (see below) 
that Hc{0'^) < 2(2 — \/2) ks 1.172, which is less than 
Hc{Rc) ~ 1.45. Therefore the numerical curves appear 
to head toward the wrong limit for i? —0 (see Fig. [T^). 

The solution of this quandary is in the role played by 
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FIG. 1: (a) Coercive field He versus disorder strength R for 
the 3d driven RFIM at T = 0 (simulation results). Triangles: 
full PBC. Stars: full PBC with a preexisting interface. Dotted 
and dashed lines: fnll PBC with a cubic seed of side a. (b) 
Sketch of the ascending branch of the hysteresis loop in the 
magnetization-field plane; the exact (full line) and effective 
(dashed line) coercive fields of (a) are also indicated. 

rare events which induce extremely strong finite-size ef¬ 
fects. The true coercive field Hc{R) is actually mono- 
tonically increasing with R, contrary to what found in 
the above simulations. For R strictly less than Rc 2^ it 
coincides with the depinning threshold of an infinite 
surface along the most difficult depinning direction (on 
depinning, see Refs. [l3, [1^0). To see this, consider 
a field H just slightly larger than the threshold value for 
the depinning transition (henceforth we focus on the as¬ 
cending branch of the hysteresis loop: see Fig. lb). In 
this regime compact droplets of up spins of linear size 
r larger than the longitudinal correlation length for a 
pinned interface, ^ {H — ^ g^j-g unsta¬ 

ble. Their interface ballistically expands, resulting in the 
unbounded growth of the associated droplets. This can 
be shown by estimating the probability Pexp{p) that a 
droplet can expand to infinity. By decomposing it into 
independent contributions, each one of them related to 
the probability of expansion over a length much larger 
than the transverse correlation length, ^ with ( 

the roughness exponent, one finds: Pexpip) = Y[’i^=o 
where Vn = r + nK^± with K ^ 1. This product con¬ 
verges to a finite value as long as [1—p(a;)] —>■ 0 faster than 
1/x when x —>■ oo, which is certainly the case (by using 
large-deviation analysis, valid for a; ^ 1, we actually pre¬ 
dict an exponential decrease). Since the medium of down 
spins is not characterized by any long-range correlation 
besides its magnetic order, we expect that these unstable 
droplets form a dilute gas with a density p ^ where 

c is a constant of order 1. The smallest ones among them. 
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but still unstable with probability 0(1), have a linear 
size of the order of and hence are at a typical distance 
1{H) ^ Due to their ballistic expansion, they in¬ 

vade the system and induce a transition to the positively 
magnetized (stable) state in a time 

t(H) = oc (i7-i?f P)-^ exp (c'[H - , 

v[H) ^ ’ 

where v oc (H — is the velocity of an unpinned 

interface close to the transition [3 and c' is a constant 
of order one. This implies that as long a.s H > the 
negatively magnetized state is unstable. Since we do not 
expect more efficient destabilizing processes in this disor¬ 
der regime, we obtain that He = i.e., the coercive 

field coincides with the depinning threshold of an infinite 
interface. This conclusion is similar to that reached in 
Ref. from field-theoretical arguments. 

We can complete this theoretical analysis by addressing 
the issue of the potential criticality for H —>■ Hc{R)~ ■ For 
concreteness, we focus on the susceptibility x = dm/dH 
and analyze whether it diverges when approaching the 
spinodal. When changing the applied field from H to 
H + 6H below the depinning threshold, the only criti¬ 
cal behavior is due to droplets of size r > ^\\{H) that 
expand through avalanches over a scale ^±{H). From 
the known behavior near a depinning transition (^ . , 

each portion of linear size ^\\iH) of the boundary of a 
droplet should advance over a volume that in average 
scales as where r is the critical exponent 

of the avalanche-size distribution. As and ^_l diverge 
when H —>■ H~, avalanches take place on all scales. How¬ 
ever their density is exponentially suppressed since the 
droplet probability goes as exp(—cr'^) (see above). This 
leads to a singular contribution akin to that found in 
Griffiths phases : Thermodynamic quantities are sin¬ 
gular but display only essential singularities and do not 
diverge at the spinodal. For instance, the susceptibility 
should vary as 


x{H)-Xre,{H)^ [ 

OC {He - 

( 2 ) 

where w = i^depUd - 1 -b Cdep)(2 - r) -b 1], Xreg denotes 
a regular, finite, contribution; x{H) thus goes to a finite 
value when H —>■ H~ , albeit in a nonanalytical way. 

We are now in a position to explain the simulation re¬ 
sults and the role played by finite-size effects. For system 
sizes L ^ i{H) ^ a droplet of linear size at 

least ^11 is typically present in the system at the applied 
field H and the physics is controlled by its (ballistic) ex¬ 
pansion. On the other hand, for L <C i{H), no unstable 
droplets are present and the negatively magnetized state 
is stable even though H > He- Depinning is in some sense 


postponed as it requires larger applied fields. Due to the 

extremely fast rise of i{H), as soon as starts to grow, 

one observes in practice a “transition” at a field He^d{L), 

which is defined by £(i?|-^'^(L)) ~ L. Its L-dependence 

_ 1 

is only logarithmic, {L) He + c(lnL) , and 

is therefore hard to detect. 

To confirm the above interpretation, we have first stud¬ 
ied the RFIM at T = 0 on a cubic lattice of size L by 
computer simulatio n ( with 20 < L < 60 and full PBC; 
details in the SM [39j), but now starting from an ini¬ 
tial condition with a plane of up spins when all other 
spins in the bulk are down. This amounts to studying 
the depinning of an infinite interface [^ . The resulting 
coercive field He{R) corresponds to the bottom curve in 
Fig. Hh. The presence of a preexisting interface lowers 
the value of the coercive field. The latter is now found to 
decrease with decreasing disorder strength. (Note that, 
except in the vicinity of the critical endpoint. Re ~ 2.2, 
where bulk and interface mechanisms also become more 
and more intertwined, the results are weakly dependent 
on system size for L > 30.) On the basis of our previous 
arguments, the He line that we now find is thus the cor¬ 
rect spinodal transition line. 

To further investigate the role played by droplets in in¬ 
ducing finite-size effects, we have run simulations of the 
model with full PBC in the presence of an initial seed (a 
droplet) in the form of a cube of up spins in the bulk of 
down spins. As with the preexisting interface, by doing 
this we put by hand the rare events needed to induce the 
transition. We have then determined the effective coer¬ 
cive field He^^ (cr, R) at which a macroscopic avalanche 
first takes place as a function of the cube side a and of 
the disorder R (for L = 60, see SM [s^). The results 
are shown in Fig. [I] The curves interpolate between the 
full PBC in the absence of seed and the transition line 
with a preexisting interface. They can be considered as 
“iso-cube-size” crossover lines: They would be observed 
in samples of linear size L containing droplets no larger 
than (T. On the basis of our above arguments, the sharp 
crossover that resembles a transition then takes place 
when ^11 ~ cr, which leads to: He^d{(j) ~ TJe + 

[for a L £{H^f^)]. We have checked that this law 
is indeed compatible with our numerical data: see the 
SM [ 3 ^. Our numerical results further indicate that the 
typical values of L reachable in simulations {L ^ 100) 
correspond to droplet sizes of a few units cr < 3, as we 
also find by direct inspection. 


Having clarified the critical behavior and the finite-size 
effects in d = 3, we turn to the issue of the potential criti¬ 
cality in high dimension and the relation with MF results. 
On the fully connected lattice and on Bethe lattices 311 
the susceptibility x = dm/dH diverges when approach¬ 
ing the spinodal as [He{R) — H]~^. This is associated 
with a scale-free, power-law, distribution of avalanches 
(with the same exponent r = 3/2 as at the critical end- 
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FIG. 2: Coercive field Hc/{Jd) versus disorder strength R 
for the driven RFIM at T = 0 in a large but finite d (a) 
and for d = oo (b). (a): Long dashed: Crossover associated 
with small droplets having cr ~ 1. Full: Actual transition 
associated with depinning and rare large droplets, (b) Full: 
Actual transition induced by purely local droplets. Compare 
also with Fig. la. 


point Rc 0)- Perturbative arguments suggest that this 
behavior persists as long as d > due = 8, but our results 
which hold in any d invalidate this conclusion. 

What is the effect of the rare events discussed above in 
high dimension? Considering the i? —>■ O'*" limit is partic¬ 
ularly instructive. The instability due to the depinning 
and the subsequent expansion of rare droplets implies 
that the coercive field in a hyper-cubic lattice in d dimen¬ 
sions is less than or equal to the threshold for the depin¬ 
ning of a (d— l)-manifold (along the most difficult direc¬ 
tion) in d dimensions: 

Moreover, a sufficient condition for depinning is that an 
infinite hyper-plane of -1-1 spins in a bulk of —1 spins 
is able to move by one lattice spacing. Consider now 
the adjacent hyper-plane that the infinite manifold has 
to invade. The spins belonging to it exactly form a 
(d — l)-dimensional RFIM. If the applied field is larger 
than the coercive value of this model, 
then the configuration flips to the positively magne¬ 
tized state and the manifold indeed moves. This implies 
jj{d,dep) and, in consequence, 
jj{d,RFiM)^Q+^ < The conclusion is 

that whatever d, = 

2(2 — '/2) Ri 1.172 [^. On the other hand, the anal- 

S ’; of the RFIM when d = oo (fully connected lattice) 
leads to a coercive field that increases as R decreases 
and goes to a value oc d —>■ oo when i? —>• 0+! The out¬ 
come is sketched both for a large d and for d = oo in Fig. 
[5] after having rescaled the coupling J by I /d. 

Clearly MF theory and the limit of infinite dimensions 
get the critical behavior and even the transition line com¬ 
pletely wrong. The droplet-induced avalanches are non- 
perturbative in 1 /d and are not taken into account within 
MF treatments, which instead capture only avalanches 
triggered by a single spin flip. (Think for example of a 


Bethe lattice where flipping spins on one branch does 
not influence the propagation of an avalanche on an¬ 
other branch.) From the perspective discussed above in 
terms of cubic seeds, the MF description corresponds to 
CT = I: the predicted critical line therefore corresponds 
to a crossover, the sharper the larger d and whose sig¬ 
nature can be seen in simulations with limited system 
sizes as for d = 9 in Ref. 0 , but not to a transi¬ 
tion. The actual transition is located in a very differ¬ 
ent region of the parameter space. In the intermediate 
droplet-dominated regime (see Fig. [Da) both the time- 
scale to evolve from one magnetic state to the other and 
the distance between droplets diverge when d —>■ oo: This 
regime is non-perturbative in 1 /d. (Note that there is no 
sharp change of the crossover in d = 8, but that the per¬ 
turbative fluctuations may nonetheless have a stronger 
influence below this dimension.) 


What are the consequences of our findings on the cur¬ 
rent understanding of glassy dynamics? On the basis 
of the mapping that relates the dynamical transition of 
glassy systems to the spinodal of the RFIM ll|, we can 
conclude that the physics described by the MCT is one 
of avalanches of motion created by a localized rearrange¬ 
ment. (The analog for glassy dynamics of a spin flip 
in the RFIM is a localized relaxation event.) From this 
point of view there is a clear connection with the phe¬ 
nomenon of dynamical facilitation (33 - [36j| . With this 
in mind, it would be worth directly establishing from 
simulations and experiments whether facilitation and 
avalanches of motion are also generated, before thermally 
activated processes take over, by rare droplets whose size 
grow as one lowers the temperature. Finally we stress 
that our results imply that even when activated events 
are neglected the MCT singularity can only have the 
meaning of a crossover, which is sharper the higher d. 
This has also been recently argued by Rizzo , but con¬ 
trary to his conclusion we find that disorder does not de¬ 
stroy the spinodal transition altogether; yet, it drastically 
changes its nature and its location via nonperturbative 
phenomena. The MCT description then becomes invalid 
whenever thermal activation or droplets with ct ^ I take 
over 0. 

We conclude by noting that the instability or failure of a 
metastable state in the presence of quenched disorder is 
a quite general phenomenon with applications from ma¬ 
terials to social science. The implication of our results 
on the role of rare relaxation events for these cases is 
certainly an issue worth being addressed in the future. 
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Supplemental Material: Spinodals with Disorder: from Avalanches in Random Magnets to Glassy Dynamics 


In this Supplemental Material we first discuss the ef¬ 
fect of the chosen zero-temperature dynamics, then give 
some details on the simulation, and finally expand on 
the results that we have obtained with the full boundary 
conditions and by introducing cubic seeds in the initial 
state of the system. 


Zero-temperature dynamics 

The quasi-statically driven RFIM at zero temperature 
has been studied with two different dynamics. One was 
introduced by Sethna and coworkers 0,1 and corre¬ 
sponds to the zero temperature limit of the Metropo¬ 
lis dynamics. Robbins and coworkers 03 introduced 
another one: a domain-wall motion corresponding to a 
domain-growth dynamics, in which the spins may flip 
under the same condition as discussed in the main text, 
but with the restriction that their location is at the in¬ 
terface between the domains of up spins and down spins; 
in this case, a preexisting interface is introduced at the 
initial time. 

In the domain-growth dynamics, the critical field is 
that of the depinning of the interface. For the 3 —d RFIM 
with gaussian random fields studied in the present work, 
Koiler and Robbins 3 found that below a (tri) critical 
value Rtc ~ 2.52, the growing domain is compact (the 
magnetization m therefore jumps discontinuously at the 
transition) and the interface is self-afHne whereas above 
Rtc one finds a self-similar percolation-like mechanism 
that results in a continuous evolution of m{H). The re¬ 
sults are displayed together with those of the present 
study in Fig. [31 Note that the critical/coercive fields 
for the two dynamics seem to coincide within numeri¬ 
cal accuracy below a disorder strength R Ri 1.5. Above 
this value, the Sethna et al. dynamics allows nucleation 
events in the bulk and the preexisting interface therefore 
moves in a more favorable environment, with a reduced 
effective disorder. In consequence, depinning takes place 
at a lower value of H than in the domain-growth dynam¬ 
ics. As one approaches the critical endpoint Rc ~ 2.2, 
bulk and interface mechanisms become more and more 
intertwined. As already pointed out 0ii , the crit¬ 
ical behavior near Rc is then different from that of the 
domain-growth dynamics around Rtc- 


Simulation details 

In our simulation, to accelerate the dynamics we have 
used a faster-than-the-clock algorithm y. We first look 
at the spins that are unstable. We use a rejection-free 
algorithm in which we sample the random time interval 



FIG. 3: (a) Coercive field He versus disorder strength R for 
the 3-d driven RFIM at T = 0 (simulation results). In addi¬ 
tion to our data shown in Fig. la of the main text we have 
plotted as a dotted line the depinning transition line obtained 
with the domain-growth dynamics in Ref. 0]. Long dashed 
line: our data with full PBC (L=90). Full line: our data with 
full PBC and a preexisting interface. 


for the next spin to be flipped. The waiting time for flip¬ 
ping an unstable spin is a Poisson variable whose mean 
is determined by the fraction of unstable spins. We then 
choose a spin at random from the unstable ones and flip it 
with probability one. We next re-compute the local fields 
on the neighboring spins and update the list of unstable 
spins. Since this is a local process, the computation cost 
is small and the algorithm is much faster than with the 
standard Metropolis one. We have checked that the out¬ 
put obtained in this way is the same as that obtained 
with the plain Metropolis algorithm. 

We have investigated three types of geometry and ini¬ 
tial conditions for the 3-d RFIM on a L x L x L cubic 
lattice with periodic boundary conditions (PBC) in all 
directions: (I) starting with all spins down, (2) starting 
with a preexisting plane of up spins at half-height in the 
sample and a bulk of down spins, and (3) starting with 
a cubic seed of up spins in a bulk of down spins. 

For case (I), we have considered random-field strength 
R from 1 to around 2.2, the latter corresponding to the 
critical point which we have checked from finite-size scal¬ 
ing analysis. We have taken L from 20 to 120 and a num¬ 
ber of samples J\f (i.e., realizations of the random fields) 
depending on the value of i?: For R = 1, J\f = 2400 
for L = 80, 90,100 and M = 6000 for smaller values of 
L. For 1 < i? < 1.5 we used L = 90 and M = 2400. 
For larger values of i? a smaller number of samples, typ- 
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ically 1200 was found sufficient. Note that even with the 
faster-than-the-clock algorithm, a simulation for R = 1 
and L = 120 takes of the order of two days on a 2.6 GHz 
computer. 

For the case (2) we have taken L from 20 to 60 with 
Af = 2400 for the smaller systems and 1200 for the larger 
ones. 

Finally, for case (3), we first checked that finite-size 
effects due to a too small ratio between the system size 
and the seed size are negligible when Lja < 3. We have 
then studied a from 3 to 20 with L = 60 and N = 1200. 

Note that the intrinsic simulation errors are very small 
and always less than the symbols used in the figures 

Data Analysis 



FIG. 4: Effective coercive field He for the 3-ci driven RFIM 
at T = 0 with full PBC and no preexisting plane or seed, (a) 
He versus disorder strength R. The triangles correspond to 
L = 90 and the circles to the result of a pseudo finite-size 
scaling analysis illustrated in (b). (b) Fit of Hc,l for 7? = 1 
and L from 20 to 120 to the form Hc,l = + AL~^ with 

He = 2.11, A = 2.78, and x = 0.47. 


cube side a for a large system size L = 60 and sev¬ 
eral values of the disorder strength R. The difference 
He^f{a, L, R) — (R) is plotted versus cr on a log-log 

scale in Fig. [5l For \na < 1.6 the effective coercive field 
is essentially constant and is given by the value in the 
absence of seed and for a ^ L one has severe finite-size 
effects. However, one observes in between a power law. 
The effective exponent is smaller but roughly compatible, 
considering the uncertainty of the present determination, 
with Xjvdep in d = 3: we find Ijvdep ~ 1-05 whereas the 
expected value is l/vdep ~ 1-14 — 1.3.0 



In a 

FIG. 5: T = 0 driven RFIM in d = 3 with a cubic seed of side 
a. The system size is L = 60. Difference with the coercive 
field in the presence of a preexisting interface as a function 
of CT on a log-log plot. The slope of the curves is about 1.05, 
within 20% of the estimated value of l/zz^ep ~ 1-14 — 1.3 0. 


When the system is studied with full boundary condi¬ 
tions [case (1)], even with system sizes as large as L = 120 
one observes an effective coercive field as exponentially 
rare droplets cannot be found in practice. In Fig. we 
show Hc{R) obtained for L = 90 as well as the result of 
a pseudo finite-size scaling analysis in which we fit the 
effective coercive field for each value of the random-field 
strength R to the form Hc.l = + AL~^. As seen 

from Fig. |3 |d the fit is good over the range of L studied 
(from 20 to 120) and gives values of comparable to 
those obtained for the largest values of L. The pseudo 
finite-size scaling procedure should however break down 
for (much) larger values of L since is not actually 

constant but rather decreases logarithmically with L, as 
discussed in the main text. 

As explained in the main text, in the case (3) where 
a cubic seed is present, we have studied the variation 
of the effective coercive field {cr, L, R) with the 


* Electronic address: sarojnandi@gmail.com 
t Electronic address: giulio.biroli@cea.fr 

* Electronic address: tarjus@lptmc.jussieu.fr 

[1] J. P. Sethna et ai, Phys. Rev. Lett. 70, 3347 (1993). 

[2] O. Perkovic, K. Dahmen and J. P. Sethna, Phys. Rev. B 
59, 6106 (1999). 

[3] H. Ji and M. O. Robbins, Phys. Rev. A 44, 2538 (1991); 
Phys. Rev. B 46 ,14519 (1992). 

[4] B. Koiler and M. O. Robbins, Phys. Rev. B 62, 5771 

( 2000 ). 

[5] K. Dahmen and J. P. Sethna, Phys. Rev. Lett. 71, 3222 
(1993); Phys. Rev. B 53, 14872 (1996). 

[6] W. Krauth, Statistical Mechanics: Algorithms and Com¬ 
putations (Oxford University Press, Oxford,2006). 

[7] P. Le Doussal, K. J. Wiese, and P. Chauve, Phys. Rev. B 
66, 174201 (2002); Phys. Rev. E 69, 026112 (2004). 













